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Abstract 

We propose new jump-adapted weak approximation schemes for stochas- 
tic differential equations driven by pure-jump Levy processes. The idea is 
to replace the driving Levy process Z with a finite intensity process which 
has the same Levy measure outside a neighborhood of zero and matches 
a given number of moments of Z. By matching 3 moments we construct a 
scheme which works for all Levy measures and is superior to the existing 
approaches both in terms of convergence rates and easiness of implemen- 
tation. In the case of Levy processes with stable-like behavior of small 
jumps, we construct schemes with arbitrarily high rates of convergence by 
matching a sufficiently large number of moments. 
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1 Introduction 

Let Z be a d-dimensional Levy process without diffusion component, that is, 
Zt = jt+ f I yN{dy,ds)+ I / yN{dy,ds), t€[0,l]. 

Jo J\y\<l Jo J\y\>l 

Here 7 e M'', N is a Poisson random measure on R'' x [0, 00) with intensity v 
satisfying / lA\\y\\'^h'{dy) < 00 and N{dy,ds) = N{dy,ds) — v{dy)ds denotes the 
compensated version of N. We study the case when i/(M'') = 00, that is, there is 
an infinite number of jumps in every interval of nonzero length a.s. Further, let 
X be an M"-valued adapted stochastic process, unique solution of the stochastic 
differential equation 

Xt = Xo+ [ h{X,_)dZs, t€[0,l], (1) 
Jo 
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where h is an m x d matrix. 

In this article we are interested in the numerical evaluation of E[f{Xi)] for a 
sufficiently smooth function / by Monte Carlo, via discretization and simulation 
of the process X. We propose new weak approximation algorithms for ([T]) and 
study their rate of convergence. 

The traditional method to simulate X is to use the Euler scheme with con- 
stant time step 

This method has the convergence rate [HI IS] 

\E{![X^)\~E{j[^Xl)\\<- 

n 

but suffers from two difficulties: first, for a general Levy measure there is no 
available algorithm to simulate the increments of the driving Levy process and 
second, a large jump of Z occurring between two discretization points can lead 
to an important discretization error. 

A natural idea due to Rubenthaler [11] (in the context of finite-intensity 
jump processes, this idea appears also in [H |H]), is to approximate Z with a 
compound Poisson process by replacing the small jumps with their expectation 

Zl:=-f,t+f f yN{dy,ds), 7, = 7 - / yiyidy), 

Jo J\y\>e "'e<|y|<l 

and then place discretization dates at all jump times of Z^ . 

The computational complexity of simulating a single trajectory using this 
method becomes a random variable, but the convergence rate may be com- 
puted in terms of the expected number of discretization dates, proportional to 
Ag = J^y^y^v{dy). When the jumps of Z are highly concentrated around zero, 
however, this approximation is too rough and the convergence rates can be 
arbitrarily slow. 

In [7], the authors proposed a scheme which builds on Rubenthaler's idea of 
using the times of large jumps of Z as discretization dates but achieves better 
convergence rates. Their idea is, first, to approximate the small jumps of Z with 
a suitably chosen Brownian motion, in order to match not only the first but also 
the second moment of Z, and second, to construct an approximation to the so- 
lution of the continuous SDE between the times of large jumps. Similar ideas of 
Gaussian correction were recently used in [S] in the context of multilevel Monte 
Carlo methods for the problem Q. However, although diffusion approximation 
of small jumps improves the convergence rate, there are limits on how well the 
small jumps of a Levy process can be approximated by a Brownian motion. In 
particular, the Brownian motion is a symmetric process, while a Levy process 
may be asymmetric. 

In this paper we develop new jump-adapted discretization schemes based 
on approximating the Levy process Z with a finite intensity Levy process Z"^ 
without diffusion part. Contrary to previous works, instead of simply truncat- 
ing jumps smaller than e, we construct efficient finite intensity approximations 
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which match a given number of moments of Z. These approximations are su- 
perior to the existing approaches in two ways. First, given that is a finite 
intensity Levy process, the solution to ([T]) with Z replaced by is easy to 
compute, either explicitly or with a fast numerical method, making it straight- 
forward to implement the scheme. Second, by choosing the parameters of Z^ in 
a suitable manner, one can, in principle, match an arbitrary number of moments 
of Z and obtain a discretization scheme with an arbitrarily high convergence 
rate. 

The paper is structured as follows. In Section [2j we present the main idea 
of moment matching approximations and provide a basic error bound for such 
schemes. In Section|3j we introduce our first scheme which is based on matching 
3 moments of Z and can be used for general Levy processes. For Levy processes 
with stable-like behavior of small jumps near zero, the scheme is shown to be 
rate-optimal. Finally, Section |4] shows how schemes of arbitrary order can be 
constructed by matching additional moments, once again, in the context of Levy 
processes with stable-like behavior of small jumps. 



2 Moment matching compound Poisson approx- 
imations 

Let Z'^ be a finite intensity Levy process without diffusion part approximating 
Z in a. certain sense to be defined later: 

:=7.t+ / / yN'{dy,ds), (2) 

Jo JR'' 

where N"^ is a Poisson random measure with intensity measure dt x such that 
:= ly^W^) < oo. 

In this paper we propose to approximate the process ([ij by the solution to 

dXt ^ h{Xt-)dZl Xo^Xo, (3) 

which can be computed by applying the Euler scheme at the jump times of Z'^ 
and solving the deterministic ODE dXt = h{Xt)ledt explicitly (or by a Runge- 
Kutta methocQ between these jump times. The following proposition provides 
a basic estimate for the weak error of such an approximation scheme. We impose 
the following alternative regularity assumptions on the functions / and h: 

(H„) / e C", h e C" /C^) and h^^^ are bounded for 1 < fc < ti and / z^-^v{dz) < 
oo. 

(H'„) / e C", h e C", /iC^) are bounded for 1 < /c < n, J^*^) have at most 
polynomial growth for 1 < fc < n and / |^|'^z^((iz) < oo for all fc > 1. 

^In this paper, to simplify the treatment, we assume that the ODE is solved explicitly. 
Upper bounds on the additional error introduced by the Runge-Kutta method are given in 
[3 Proposition 7]. These bounds can be made arbitrarily small by taking a Runge-Kutta 
algorithm of sufficiently high order. 
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Proposition 1. Let Z and Z be Levy processes with characteristic triplets 
(0, v, 7) and (0, i), 7) respectively, and let X and X be the corresponding so- 
lutions of SDE ([1]). Assume ^ — ^, u — v on {\\x\\ > 1}, either (Hn) or (H^) 
for n > 3 and 



Xi^ . . .Xif^y{dx) = I Xi^ . . .Xif^v{dx), 2<k<n-l, I < ik < d. (4) 



Then 

\E[f{X,)~ fiX,)]\<C f \\xr\du-di^l 

where the constant C may depend on f , g, x and v but not on v. 

Proof. To simplify notation, we give the proof in tlie case m ~ d — 1. Let 
u{t,x) = By Lemma 13 in [7 , u e Ci'"([0, 1] x M) and satisfies 

(i, x) + 7— (i, a;)/i(x) + / {u{t,x + h{x)y) — u{t,x))v {dy) 
J\v\>i 

Ou \ 

u{t, X + h{x)y) ~ u{t, (tj x)h{x)y j v{dy) = 0, (5) 



dt 



dx 



'|y|<i 
■"(1,2;) = f{x). 

Applying Ito formula under the integral sign and using ^ and Lemma 11 in 
[7] (bounds on moments of Xt) yields 

E[f[X^) ~ /(Xi)] = E[u(l, Xi) - u(0, Xo)] 

u{t, Xt + h(Xt)z) - u{t, Xt) - h{Xt)z— \ {di^e - dv)dt 



= E 
+ E 
= E 



dx 



^^u{t,Xt- + h{Xt^)z) - u{t,Xt-)^ N{dz,dt) 
' - di^jdt + remainder 



k=2 

E' [remainder], 



where in the last line we used the moment matching condition (17) and the 
remainder coming from the Taylor formula can be estimated as 



I remainder I < 



sup 

JmO<s<1 



< C sup 

0<s<l 



ax" 



\h{Xt)nzndiy, - d,^\dt 



sup / \z\^\du, - di^\ 

0<s<l 



From the Lipschitz property of h and Lemma 13 in [7], 



sup 

0<s<l 



ax" 



sup \h{X,T<C{l+ sup iXtin 

0<s<l 0<t<l 
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for some C < cxi, where p = n under (Hn) and p > n under (H^). Following 
the arguments in the proof of Lemma 11 in [7j, we get 



E[ sup \Xt\P] < C(l + exp 

0<t<l 



p/2\ 

c| |7|P+ / \zfu{dz)+[ I z^v{dz) 



for different constants C and c, where 



7 = 7+ / zi){dz) = 7 + / zv{dz) 
J\z\>l J\z\>i 

by our assumptions. Since ^^ z'^i'{dz) = J^z'^h'{dz) by assumption, and 

\z\PHdz) < [ \z\Piy{dz)+ [ \z\^iy{dz) = / \z\P:^{dz)+ [ \z\^iy{dz), 

J\Z\>1 "'|2|<1 "'|2|>1 "'|Z|<1 

it is clear that £'[supQ<f<]^ l^tl^] 1^ C for some constant C which does not 
depend on v. 

□ 



3 The 3-moment scheme 

Our first scheme is based on matching the first 3 moments of the process Z. 
Let S'^^'^ be the unit sphere in the d-dimensional space, and v{dr x d9) be a 
Levy measure on written in spherical coordinates r S [0, oo) and 6 e S"^^^ 
and satisfying Jj^ c»)xS'^-i ''^^^(.dr, dO) < oo. Denote by u the reflection of ly with 
respect to the origin defined by P{B) = i^{{x : —x € B}). We introduce two 
measures on S'^~^: 



1 

2 7|r|<e £^ 



X(de) = - / {v{dr, dO) - v{dr, dO)) 



\{de) = - / — {v{dr, dO) + v{dr, dO)) . 

2 J\r\<e £ 

The 3-moment scheme is defined by 

v^{dr,de) = v{dr,d6)lr>e + 5e{dr){\{d6) + X{de)} (6) 



le = l- / r0iy,idr,de), (7) 



where 6^ denotes a point mass at e. 
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Proposition 2 (Multidimensional 3- moment scheme). For every e > Q, is 
a finite positive measure satisfying 



XiXjv{dx) — I XiXjVe{dx) (8) 
XiXjXki'{dx) = I XiXjXkVf;{dx), i<i,j,k<d (9) 



Ae := / i^,{dx)= I i^(da;)+£"M ||a:||^j/(da;) (10) 

II >£ ■/ II a;|| <s 



[ \\x\\*\diy - dl^el < [ 
JR-i J\\, 



\\x\\ v{dx ) +eM \\x\Yv{dx), (11) 

a;||<£ ^||a:||<e 



where the last inequality is an equality if u{{x : ||.t|| = e}) = 0. 

Proof. The positivity of being straightforward, let us check ([8|. Let {ei\f^^ 
be the coordinate vectors. Then, 

XiXjV^{dx) — I r'^{9,ei){9,ej)i's{dr,d6) 

' J[0,oo)xS<*-i 

r^{9,ei}{e,ej)i^{dr,de)+ [ e^{0,e^}{e,ej}{\{de) + X{d0)} 
r^{9,ei){e,ej)i^{dr,de) + [ {0 , e,) {9 , e^) X{de) 



r {9,ei){0,ej)v{dr,d9) = I XiXjv[dx). 

(0,oo)xS<*-i JK<* 

The other equations can be checked in a similar manner. □ 
Corollary 1. Let d = 1. Then the 3-moment scheme can be written as 

Vs{dx) — v{dx)l\x\>€ + \+6^{dx) + \-5-e{dx) 

" f ] 

;v{dx) ± / -^v{dx) > 



Corollary 2 (Worst-case convergence rate). Assume (H4) or (H4). Then the 
solution X of ^ with the characteristics of given by ^ -([T]) satisfies 

\E[f{X,)-f{X,)]\=o{X-'). 

as e ^ Q. 

Proof. By Proposition [l] we need to show that 

limAe / \\x\\'^\dv ~ dv,\ =0. 
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By Proposition [2j 

limAe / ||a;||^|(ii^ — (i^gj 




<lim< v{dx)+e-' \\xYv{dx) 

■^^0 I J\\x\\>e "'l|2:||<e 



x\\ <e J \\x\\ <£ 




<21ime2<' / v{dx) + e-'' I \\x\\^v{dx))l \\xfv{dx) 

21inie2 / v{dx) I \\x\\^v{dx) 
<2 I \\xfv{dx)Vmie^ v{dx)=0, 

Js.'i J\\x\\>e 

where in the last line the dominated convergence theorem was used. □ 

In many parametric or semiparametric models, the Levy measure has a sin- 
gularity of type near zero. This is the case for stable processes, tempered 
stable processes [TU], normal inverse Gaussian process [T], CGMY [3] and other 
models. Stable-like behavior of small jumps is a standard assumption for the 
analysis of asymptotic behavior of Levy processes in many contexts, and in our 
problem as well, this property allows to obtain a more precise estimate of the 
convergence rate. We shall impose the following assumption, which does not 
require the Levy measure to have a density: 

(H - a) There exist C > and a e (0, 2) such thai[^ 

l{r) - Cr'" as r ^ (12) 

where l{r) := J^^^^^^^,iy{dx). 

Corollary 3 (Stable-like behavior). Assume (H — a) and (H4) or (H4). Then 
the solution X of ^ with the characteristics of given by Q -Q satisfies 



\E[f{X,)^fiX,)]\=o[\l-^) 



Proof. Under (H ~ a), by integration parts we get that for all n > 2, 

x|l"z/(dT) - -^r""" asr^O. 
x\\<r n-a 

Therefore, under this assumption, 

2C f i,ii4|,., , , I ^ Ca(6-2a) 



Ae e-" and / \\x\ndiy - di^J < 1 ^7 as e ^ 0, 

2-a Vr. ~ (2-a)(4-a) 

from which the result follows directly. □ 



■^Throughout this paper we write / ~ <; if hm = 1 and / < g if limsup < 1. 
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Rate-optimality of the 3-moment scheme From Proposition [T] we know 
that under the assumption (H4) or (H4), the approximation error of a scheme 
of the form can be measured in terms of the 4-th absolute moment of 

the difference of Levy measures. We introduce the class of Levy measures on 
M'' with intensity bounded by N: 

= {v Levy measure on M^, < N}. 

The class of Levy measures with intensity bounded by is then denoted by 
Ai^^ , and the smallest possible error achieved by any measure within this class 
is bounded from below by a constant times inf^/g^Ae Jj^^ ||a;||'^|di' — dj^'\. The 
next result shows that as £ — ?> 0, the error achieved by the 3-moment scheme ly^ 
differs from this lower bound by at most a constant multiplicative factor . 

Proposition 3. Assume (H — a) and let he given by ^. Then, 

Proof. Step 1. Let us first compute 

£n:= inf / \\x\\^\dv ~ du' \. (13) 

For v' ^ ^\ci v' ^ v'^ + v'^ where v'^ is absolutely continuous with respect to 
V and v'g is sin gular. Then v'^{M.'^) < N and 

\\x\\^\diy-diy'\= f \\xf\diy-diy'^^\+ f \\x\\Uiy'^. 



Therefore, the minimization in ( 13 1 can be restricted to measures v' which are 



absolutely continuous with respect to ly, or, in other words, 
Ej^ = inf / \\x\\*\l - X{x)\iy{dx), 



where the inf is taken over all measurable functions A : R'' M"*" such that 
/gd \{x)v{dx) < N. By a similar argument, one can show that it is sufficient 
to consider only functions A : M*^ — >■ [0, 1]. Given such a function X{r,9), the 
spherically symmetric function 



A(r) 



Jgd-1 ^{dr, dO) 



leads to the same values of the intensity and the minimization functional. There- 
fore, letting v{dr) := J^d-i v{dr,d9), 

/>00 POO 

£n^ inf / r^{l- X{r))i}{dr) s. t. / X{r)i}{dr) < N. (14) 

0<A<lJo Jo 
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For every e > 0, 

En> inf (/ r^{l-\{r))v{dr) + e^ ( [ \{r)v{dr) ~ N 

0<A<1 LJo \Jo 

The inf in the right-hand side can be computed pointwise and is attained by 
A(r) — lr>e + f^^r=e for any fi G [0, 1]. Let e{N) and At(iV) be such that 

z>((e(7V), oo)) + ^liN)^>{{e{N)}) = N. 

Such a e(-/V) can always be determined uniquely and /i(-/V) is determined uniquely 
if iy{{e{N)}) > 0. It follows that A(r) = lr>e(7V) + /^(-^)lr=e(Af) is a minimizer 
for ( 14 1 and therefore 

£n= f M^Hdx) + (1 - M)eV({x : = e}), 

^ 1 1 x 1 1 < e 

where e and /x are solutions of 

: ||a;|| > e}) + : ||a;|| = e}) = N. 

Step 2. For every e > 0, let e(e) and ^(e) be solutions of 

v{{x : \\x\\ > e(e)}) + fi{e)i^i{x : \\x\\ = e(e)}) = A^. 

It is clear that e(e) — > as e — > and after some straightforward computations 
using the assumption (H — a) we get that 

e(e) 
lim — ' 

e-H-O e 



Then, 



e(e)4-" 



X hm ■ 



i!|r.i|<e(e) + " K^M^Vl^ii^ ■ M = e(£)}) 

Under (H — a) the three limits are easily computed and we finally get 

j.,My.-...i /^y/' (15) 

□ 

Remark 1 . The constant (3 — a) I J > 1 appearing in the right-hand side 

of (151 cannot be interpreted as a "measure of suboptimality" of the 3-moment 
scheme, but only as a rough upper bound, because in the optimization problem 
( 13 ) the moment- matching constraints were not imposed (if they were, it would 
not be possible to solve the problem explicitly). On the other hand, the fact that 
this constant is unbounded as a — ?> 2 suggests that such a rate-optimality result 
cannot be shown for general Levy measures without imposing the assumption 
(H-a). 
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Numerical illustration We shall now illustrate the theoretical results on a 
concrete example of a SDE driven by a normal inverse Gaussian (NIG) process 
[1], whose characteristic function is 

0f (u) := £;[e™^'] = exp (V"^-(/3-m)2 - ^ a'^ - /32^ | , 

where a > 0, /3 e (— ct, a) and (5 > are parameters. The Levy density is given 

by 

v{x) = — , 

where K is the modified Bessel function of the second kind. The NIG process has 
stable-like behavior of small jumps with v{x) ~ ^^yr, a; — >■ (which means that 
(H — a) is satisfied with a = 1), and exponential tail decay. The increments of 
the NIG process can be simulated explicitly (see [4j algorithms 6.9 and 6.10]), 
which enables us to compare our jump-adapted algorithm with the classical 
Euler scheme. 

For the numerical example we solve the one-dimensional SDE 

dXt — sin(aXf )c?Zi, 

where Z is the NIG Levy process (with drift adjusted to have E\Z^ — 0). The 
solution of the corresponding deterministic ODE 

dXt = sin(oX()di, = x 

is given explicitly by 

^, , 1 1 + cos(aa;) — e2°*(l — cos(ax)) 

Xt — &{t; x) = - arccos — r 5-— — 

a 1 -I- cos(ax) -I- e^"'(l — cos(ax)) 

Figure [T] presents the approximation errors for evaluating the functional 
E[{Xi — 1)+] by Monte-Carlo using the 3-moment scheme described in this 
section (marked with crosses), the diffusion approximation of [7] (circles) and 
the classical Euler scheme (diamonds). The parameter values are a = 0.5, 
6 = OA, K = 0.6, a = 5 and Xq = 1. For each scheme we plot the logarithm of 
the approximation error as function of the logarithm of the computational cost 
(time needed to simulate 10^ trajectories). The curves are obtained by varying 
the truncation parameter e for the two jump-adapted schemes and by varying 
the discretization time step for the Euler scheme. 

The approximation error for the Euler scheme is a straight line with slope 
corresponding to the theoretical convergence rate of ^. The graph for the 3- 
moment scheme seems to confirm the theoretical convergence rate of Xj^; the 
scheme is much faster than the other two and the corresponding curve quickly 
drops below the dotted line which symbolizes the level of the statistical error. 
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Figure 1: Approximation errors for the 3- moment scheme (crosses), the scheme 
of [7] (circles) and the Euler scheme (diamonds). The horizontal dotted line 
corresponds to the logarithm of the two standard deviations of the Monte Carlo 
estimator (the standard deviation is about the same for the three schemes and 
independent of the discretization step) ; everything that is below the dotted line 
is Monte Carlo noise. 

4 High order schemes for stable-Uke Levy pro- 
cesses 

In this section, we develop schemes of arbitrary order for Levy processes with 
stable-like behavior of small jumps. Throughout this section, we take d = 1 
and let Z be a Levy process with characteristic triplet (0, v, 7) satisfying the 
following refined version of (H — a): 

(H' - a) There exist, c+ > 0, c_ > with c+ + c_ > and a £ (0, 2) such that 

poo I* — 6 

/ v{dx) ^ c+e^" and / v{dx) ^ c_e"" as £ 4- 

J e J —00 

Introduce the probability measure 

X (2 - a)|a;|i""(c+lo<x<i + c_l_i<j,<o) 
M (x) — p-^ — . (16) 

Let n > and e > 0. The high-order scheme for the stochastic differential 
equation ([T]) based on n + 2 moments and truncation level e is constructed as 
follows: 
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1. Find a discrete probability measure /i = Y^^^q o-i^xi with 

x^ix{dx)^ I x^li*{dx), l<k<n, (17) 



such that xo < xi < ■ ■ ■ < Xn, Xi ^ Q for all i and a* > for all i. 

2. Compute the coefficients {af} by solving the linear system 

<7g af a;*'e'' = / x^+^u{dx), k = 0,...,n, cr^ = / x^v{dx). 

3. The high-order scheme is defined by 



/ T \ / 7 \ 2 X ^ ^7 ^SX (dX^ I / 7 \ 

iye{dx) = iy{dx)l\^\^^ + 2^2 ' 7e = 7 - / z^^eia^^) 

i=o 



/ 



1^1: 

(18) 

Remark 2. The first step in implementing the scheme is to solve the moment- 



matching problem (17 1 for measure ^* . The existence of at least one solution 
to this problem with xq < xi < . . .Xn and a* > for all i is guaranteed by 
the classical Caratheodory's theorem, but this problem admits, in general, an 
infinite number of solutions. Here we impose the additional condition 7^ 
and a* > for all i, which should be checked on a case by case basis in concrete 
realizations of the scheme (see Example [l]) . 

Remark 3. It is easy to see that the measure 

.*(dx):=(a:)^^%^, (Kr-.^f x^v*(z)dz. 



4=0 



\x\<E 



matches the moments of orders 2,...,n -I- 2 of v*(x)\x\<e^ where v* is the 
measure given by 

Q;c+la;>o + ac_l2:<o 

that is, V* satisfies the assumption (H' — a) with equalities instead of equiva- 
lences. The idea of the method is to replace the coefficients {a*} with a different 
set of coefficients while keeping the same points \xi\ to obtain a measure which 
matches the moments of v{x)\\r^\<^^. Therefore, the points \xi\ do not depend 
on the truncation parameter e while the coefficients {af } depend on it. 

Example 1 . As an example we provide a possible solution of the moment match- 
ing problem for n = 3, which leads to a 5-moment scheme (matching 3 moments 
of n* or 5 moments of the Levy process) . We assume that fi* has mass both on 
the positive and the negative half-line: c+c_ > 0. 
The moments of /j* are given by 

m.= ,^:" (p+(-l)ni-p)), p:=^^. 
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It is then convenient to look for the discrete measure matching the first 3 mo- 
ments of /i* in the form 

fl = (1 - p){p5^,, + {1 - p)5^,,) + p{{l - p)6,, +pS,,), (19) 

where p € (0, 1), < £i < £2 are parameters to be identified from the moment 
conditions 

(i~p)gj:-+p£^^ ^~" , fc = i,2,3. 

k + 2 — a 

For the purpose of solving this system of equations, let f be a random variable 
such that P[£ — 62] — p = 1 — P[S — Si]. From the moment conditions, we get: 

e:=E[£]^l^, a^:^Varg= ^\ f (20) 
6 — a (4 — q;)(3 — a)^ 

E[{£ - E[£]f] a-1 [I^ 
On the other hand, the skewness s can be directly linked to the weight p: 



" P= o - oSign(s)\/^— , (22) 



and the parameters Si and S2 can be linked to e, a and p: 



s.^e^a^l^^, e2^e + a^^. (23) 

The dependence of ei, £2 and p on a is shown in Figure [2} it is clear from 
the graph that the constraints p G (0, 1) and < £1 < £2 are satisfied for all 
a £ (0,2): therefore, equations (19 [23) define a 4-atom probability measure 
which matches the first 3 moments of /i* . 



Proposition 4. Let {xi} be fixed according to (17). There exists £0 > such 
that for all e < So, is a positive measure satisfying 

x''v{dx) = / x^v^{dx), 2<k<n + 2 (24) 

There exist positive constants Ci and C2 such that 

= i/,(M) ~ Ci£~", / |a;r+3|di/-dt^e| < C2£"+^^" as £ ^ 0. 



Corollary 4. Assume (H„+3) or (H^_,_3). Then the solution X of ^ with 
characteristics of given by ( 18 I satisfies 



\E[f{X^) - /(Xi 



OA 
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alpha 



Figure 2: Solution of the moment matching problem for 3 moments (see Example 

0. 



of Proposition The moment conditions ( 24 1 hold by construction. Using in- 
tegration by parts, we compute 

+ ass^OforA>0. 

\z\<e 2 + k-a 

Therefore, 

hm 4^ [ z-^Mdz) = ^l~''}^'^ty^''i = [ x',*idx). 



e^o t7|efc 7|^l<^ (2 + fc-a)(c++c_) 

Since the matrix Mij = {xjY, 0<i<n, 0<j<nis invertible (Vandermonde 
matrix), this impUes that hm£_j.o af = a*. Therefore, there exists Eq > such 
that for all e < Eq, a1 > for aU i and is a positive measure. 
We next compute: 



_ 2 

|x|>£ t J\x\>e t a;^ 

( n ^ \ 



2 



/ \x\^+''\dv - dv,\ < [ 

Js. J\x 



n+l 



|x|<e 

"+3-"(c++c_) 



3 + fc — a 2 — a 



E 

i=0 



□ 
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